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SYNOPSIS 


Solution the problem determining the frictional losses natural streams 
usually depends the judgment and experience the engineer. rational 
approach this problem offered this paper. Exisitng theory friction 
open channels reviewed, and logical derivation leads the development 
formulas and charts determine the relationship between stage and discharge 
given reach river. interdependent treatment the hydraulic 
characteristics flow and the sedimentary characteristics the stream 
presented. Friction losses due grain roughness, channel irregularities, and 
other factors are studied and formulated, and suggestions are made for practical 
application proposed rational method friction evaluation. 


INTRODUCTION 


dealing with natural streams with movable beds, the engineer con- 
fronted with the problem evaluating frictional losses. Unless adequate and 
complete data are available guide him, too often not the case, his evaluation 
frictional effects based almost entirely upon judgment. The element 
doubt inevitably present forces the engineer conservative. The degree 
conservatism exercised reflected project costs. the element doubt 
and confusion concerning frictional losses natural streams could reduced 
the lowest practicable limit, substantial savings the costs many hydraulic 
projects undoubtedly would result. 

Although many data concerning the roughness values natural streams have 
been accumulated and published, the fact remains that, estimating the rough- 
ness value particular stream for which roughness data are available, there 


are likely many different estimates there are estimators. Even when 
comments are invited for publication; the last discussion should submitted 
January, 1951. 
Associate Prof. Mech. Eng., Univ. California, Berkeley, Calif. 
Hydr. Engr., Omaha Dist., Corps Engrs., Army. 
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CHANNEL ROUGHNESS 


roughness data are available for particular stream, often these data are appli- 
cable only within very narrow range discharge and for particular location 
the stream. not uncommon for many engineers think natural 
river channel having single roughness value for all discharges. Even some 
those engineers who feel that the roughness value (as applied Manning’s 
Kutter’s formula) varies with the discharge are compelled, lack data 
stream. 

evident that new and more rational approach the evaluation 
frictional losses natural channels needed. The writers believe that real- 
istic appraisal frictional losses streams with movable beds possible only 
when the sedimentary characteristics the stream and the usual hydraulic 
characteristics flow are treated interdependently. With regard the sedi- 
mentary characteristics, frictional losses natural streams depend the size 
and gradation sediments composing the stream bed surface and the shape, 
size, and spacing irregularities within the boundaries flow. because 
the lack these irregularities that laboratory roughness values derived under 
hydrodynamically similar conditions often not agree with the roughness 
values natural streams. probable that these irregularities are function 
the sediment transport which certainly function flow conditions. The 
interdependence the sedimentary and hydraulic conditions flow obvious. 
the functional relationship between the irregularities the bed and the sedi- 
ment transport can expressed dimensionless form and evaluated means 
experimental data, will possible then for the practicing engineer use 
this relationship means evaluating the frictional losses natural streams 
with considerably more confidence than the past. 


THEORY 


Turbulence theories have substantiated beyond doubt that the distribution 
flow velocities real fluid (with internal friction) very closely related the 
shear stresses the fluid. Since the flow water open stream channel 
influenced friction along its solid boundaries, the bed and the banks, the 
frictional effect the flow may measured 
the shear stresses along these boundaries. 

Fig. the total shear force along the wetted 
fined the average unit shear stress, parallel 
the direction flow, within the body water 

water, cross-sectional area has weight com- 
and the slope the water surface small. For the case uniform motion 
the algebraic sum the forces acting parallel the direction flow must equal 
zero since acceleration flow takes place. Neglecting air friction the 
open surface and noting that the total shear force acts retarding direction, 
follows that 


| 
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the unit mass density, which the acceleration due gravity, 


wSA 
(2a) 
and 


Both forms Eq. are equally valuable. The dimensions the channel ap- 
pear only the ratio A/p commonly known the radius. 
the derivations this paper will prove most helpful, however, think the 
hydraulic radius the ratio the ratio the volume the flow segment 
the corresponding frictional surface area. 

Fluid mechanics teaches that hydrodynamically smooth wall the shear 
forces are transmitted entirely the tangential shearing stresses within the 
fluid, whereas hydrodynamically rough wall part the shear forces are 
transmitted unsymmetrically distributed normal stresses the individual 
roughness protrusions the wall. the latter case the degree resistance 
flow depends great extent upon the size, shape, and spacing the individual 
irregularities. This type resistance commonly called form resistance. The 
actual flow pattern the immediate vicinity rough wall basically differ- 
ent from that near smooth wall. However, distance away from the wall 
one two roughness lengths this difference disappears entirely, and the 
two flow patterns become identical. 


NATURAL STREAM BED 


rough wall discussed far wall with surface 
characterized straight alinement and practically uniform roughness. The 
pipes which roughened gluing sand grains the inside are 
good example such straight, uniformly rough surface. The individual 
grains bed sediment natural stream also provide grain roughness which 
approximately comparable that Mr. Nikuradse’s roughened pipes, and 
the resulting frictional losses due the grain roughness alone may calculated 
much the same way are the frictional losses for pipes artificial uniform 
channels immovable beds. However, the natural stream bed surface, espe- 
cially that sediment-carrying alluvial stream, basically different from that 
Mr. Nikuradse’s experimental pipes because the surface not straight. 
often irregular due sediment bars various scales. 

The foremost problem confronting engineers who deal with natural streams 
evaluate the influence the irregularities the stream bed surface upon 
the grain roughness. very likely that these surface irregularities are 
source additional friction losses over and above those caused the grain 
roughness alone. The question can these additional friction losses de- 


scribed physically. clue the answer may found from consideration 
356, VDI-Verlag, Berlin, 1932. 


rauhen Rohren,” Nikuradse, Forschungsheft 361, Berlin, 
1933. 
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CHANNEL ROUGHNESS 


the problem involving the evaluation friction losses rough pipe line which 
has irregularities such elbows, enlargements, contractions, any other type 
transitional fittings. common knowledge that, despite the very serious 
theoretical objections due the basic nonlinearity the problem, one may 
with reasonable degree accuracy determine individually the friction losses 
each element the pipe line and simply add them obtain the total loss. 
neglecting the mutual interference flow patterns within the different ele- 
ments the pipe line, the problem thus reduced that calculating the 
individual losses. The same type approach suggested this paper for the 
determination the additional friction loss due irregularities natural 
streams with moving sediment beds. 

The friction loss due grain roughness may evaluated use the 
Manning equation any other suitable open channel formula. the 
Manning equation with the average velocity 


1.486 


the roughness factor may replaced the grain roughness diameter 
means the empirical derived Strickler: 


Eqs. Eq. may rearranged into the form: 


follows from the way Eqs. were derived that Eq. the hydraulic 
radius pertaining the grain roughness. Stated another way the ratio 
the cross-sectional area (contributing the energy that converted into turbu- 
lence along the perimeter representative grain roughness K,) and the 
perimeter Call this area and the corresponding hydraulic radius R’. 


] & 


This the Manning-Strickler equation. may replaced logarithmic 
equation that better supported theory: 


This equation gives for very wide range conditions R’/K, 2,000) 
practically the same results Eq. Eq. seems very reliable de- 

zur Frage der Geschwindigkeitsformel und der Rauhigkeitszahlen Strom und 
geschlossene Leitungen,” Strickler, Mitt. No. des Amtes fiir Wasser wirtschaft, 
Bern, Switzerland, 1923. 


“Flow Movable Bed,” Einstein, Proceedings the Second Hydraulic Conference, 
Inst. Hydr. Research, State Univ. Iowa, Iowa City, 1942, 332. 
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scribing all flows hydrodynamically rough surfaces. For such flows the 
ratio the roughness diameter the theoretical thickness the laminar 
sublayer 


which the kinematic viscosity. 
smaller than Eq. must corrected. The correction may 
made introducing factor into the equation which then becomes 


which the correction function and given Fig. 2.8 This 


Correction Factor, 


0.6 


THE LOGARITHMIC FRICTION FoRMULA 


relationship derived from Nikuradse’s measurements,‘ flows over sand 
surfaces. (In some practical applications the factor may neglected since 
its effect upon the ratio V/V’; can never more than about percent.) 
Besides these usual friction losses attributable the grain roughness the 
natural stream, there are additional friction losses whose origin connected 
some manner with the presence sand bars, sand waves, and other irregulari- 
ties the channel. For the easiest understanding how these irregularities 
generate additional friction losses, imagine them uniformly sized and 
uniformly spaced protrusions the stream bed, such piles single rocks. 
Assuming the number protrusions the surface area (Fig. 1), 


Mechanics for Hydraulic Hunter Rouse, Engineering Societies Monographs, 
McGraw-Hill Book Co., Inc., New York, Y., 1938, 243, Eq. 171. 


Fig. 123. 
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each protrusion having cross-sectional area diverting the flow, the drag 
force each protrusion may expressed as: 
2 


which the drag coefficient the protrusion, depending its shape and 
size. The total drag force, the direction flow, protrusions then 
2 


The additional shear stress the boundary surface due this total 


drag force 


The total shear force Eq. may considered include two compo- 
nents: the shear force due the grain roughness and the shear force due 
the irreguiarities Similarly, the total area the cross section 
Fig. contributing energy that converted turbulence the bed, may 
broken down into two parts and that part the total area which 
pertains the grain roughness previously defined, and that part the 
total area which contributes the energy that spent friction due the 
irregularities the channel. may now written 


Since R’, follows that 


Eliminating from Eq. and Eq. and letting then 


will noted that the derivation Eq. depends the basic assumption 
that the grain roughness and the irregularities are distributed over 
the boundary surface 

The entire expression the right side Eq. function the shape, 
size, and frequency the irregularities along the surface area Now, 
thinking the irregularities (or protrusions) sand bars, most natural 
expect them function the sediment transport along the channel bed. 
This sediment transport itself function the shear stress 7’, the sediment 
grains and the representative grain diameter The effect the variables 


(17) 


for the Transportation Eed Load,” Einstein, Transactions, ASCE, Vol. 107, 


1942, 561 
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which the unit mass density the sediment air and the unit mass 
density water. Therefore, This relationship must 
determined from actual river measurements. 

Determination Additional Friction Losses Due Channel Irregularities.— 
These losses are reflected the friction velocity Eq. 16. 
absolutely uniform channel these losses are obviously equal zero, and, 
therefore, V’’; never zero real channel but varies, depending upon the 
degree the irregularities. The effect these irregularities distort the 
flow pattern. When the discharge least, the distortion the flow pattern 
greatest (as witness the meandering natural stream low flows). the 


Missouri River Near Fort Randall, S.D. 


Missouri River Near Pierre, 
Missouri River Near Omaha, Nebr. 
Elkhorn River Near Waterloo, Nebr. 
Big Sioux River Near Akron, lowa 
Platte River Near Ashland, Nebr. 
Niobrara River Near Butte, Nebr. 
Salinas River San Lucas, Calif. 
River Near Junction, Calif. 


1.68 


discharge increases (and hence the sediment transport along the bed also in- 
creases) the distortion the flow pattern becomes less and less because the 
alinement flow becomes progressively straighter. Consequently, one may 
expect that the additional friction loss (or V’’;) diminishes the discharge 
increases. Since both and become larger with increasing follows that 
the ratio also increases, but the ratio from Eq. decreases. (Note 
that increasing transport sediment means that decreases.) The 
relationship between and has been determined from actual measure- 
ments for number rivers, most which are the Missouri River Basin. 
The results obtained are indicated Fig. Except for heavily vegetated 
streams, there appears fairly well-defined relationship between 
andy’. curve for the relationship sought has been drawn Fig. 


} 
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The supporting computation for the results shown Fig. are based the 
following procedure: 


comparatively long reach river was selected within which there 
exists gaging station where the relationship between discharge and water 
surface elevation known; 

The over-all slope the reach was determined the total drop divided 
total distance; 

Surveyed cross sections the river channel spaced more less uniformly 
throughout the reach were then superimposed upon one another the gaging 
station, using the over-all slope means bringing the cross sections to- 
gether common water surface line the gaging station; 

The average river width for this common water surface line (correlated 
with discharge the gaging station) was obtained, and then the average widths 
depth intervals referred the common line were calculated. Thus, the 
average cross section for the entire reach river was determined. 

Mechanical analyis curves were constructed for all bed surface samples 
(at least samples each cross section obtained drag bucket the 
active channel) and average, composite mechanical analysis curve was then 
determined taking the average the percentage values from each individual 
curve for each chosen grain size. 

roughness. The grain size the average mechanical analysis curve which 
65% finer was taken representative value 

Hydraulic computations were made using the average cross section, the 
over-all slope and the stage-discharge relationship the gaging station. 
Thus, for every chosen discharge the cross-sectional characteristics and 
the average section were known. The hydraulic computations are based upon 
the following assumptions: 


The water surface profile essentially straight line long reach and 
remains straight with unvarying slope the discharge changes. 

The average cross section fairly long reach does not substantially 
change shape and size the discharge varies. (There some supporting 


evidence for this assumption.) 
Because and the uniform flow equations may used with negligible 


error. 
Eqs. and were then used determine R’, with all other quantities 


known these equations. 

The total was determined A/p, and the was determined 
The latter step permissible because both and are derived 
from the same wetted perimeter and equals defined. 


10. With and known, the dimensionless values and were 


1.68 D/R’S from The ratio which represents the density 


the bed sediment under water with respect the density water may 
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assumed have value between 1.65 and 1.70 for most natural sediment mix- 
tures, with 1.68 fair average. The representative grain diameter the 
sediment that transported along the stream bed was chosen that grain 
size the average mechanical analysis curve that 35% 

11. making computations for comparatively narrow streams with steep 
banks, the friction loss these banks was calculated separately friction 
rough wall, using Manning’s formula (with and known and assigned) 
determine the hydraulic radius pertaining these banks. Then the area 


wetted perimeter the banks was deducted from the total cross- 
sectional area give A’+ With total wetted perimeter 
was then computed the quotient divided the same wetted perimeter 
that was used evaluate A’. 


Description River brief description the river reaches that 
were selected for study and pertinent data are given The data fur- 
nished Table are complete enough permit computational check 
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TABLE Data 


Test 
Reac (cU PER SEC) 


AVERAGE WATER 
VELOCITY 
Approximate location sections} slope mum mum 
Mini- Maxi- Mini- Maxi- 
mum mum mum mum 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
Ft. Randall, 35.5 1.49 145,000 1.11 6.12 1,038 2,130 
Elkhorn..... Waterloo, 9.0 5.28 1.32 7.67 225 291 
Big Akron, 29.5 1.89 550 8,550 1.14 3.46 132 184 
Butte, 1.0 12.57 1,045 1,045 1.95 1.95 804 804 
Ashland, Nebr.*. 7.38 78,000 1.64 9.14 725 1,541 
Salinas...... San Lucas, 7.0 11.60 0.61 6.67 130 2,050 
Junction, 3.0 13.20 38,300 0.84 7.72 272 1,081 
Salinas...... Paso Robles, 6.0 17.20 6,000 0.14 3.00 110 595 


Subscript refers finer average mechanical analysis curve. Fairly straight 
many large sand bars. River works for navigation have materially narrowed natural bars not 
correction made. Winding stream; banks vegetated; sand bars visible; bank correction made. Very 
many split channels; winding stream; great number sand bars; range stage for large change 
vegetation banks. Fairly straight channel divided wood strips; heavily vegetated. 


least one poirt for each stream indicated Fig. except the Elkhorn and 
Big Sioux Rivers for which bank corrections were made. 
Practical Application 
2000 Relationship Between 
and adequate hy- 
draulic data are lacking for 
1600 close determination the 
stage-discharge relationship 
particular stream, the 
curve Fig. may used 
determine this relationship. 
How this may done best 
illustrated example. 
Suppose desired de- 
termine the rating curve for 
stream whose characteristics 
(having been determined 
surveys and application 
procedures previously de- 
scribed) are follows: (1) 
Stage (or Depth), Feet Slope 
cross-sectional dimensions 
the average cross section 
the channel fairly long reach function stage indicated Fig. 
(3) the grain size from the average mechanical analysis curve within the chosen 
reach which 65% finer 0.50 mm, and the grain size from the same 


1200 


Hydraulic Radius, 


800 


Area, Square Feet 


400 
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SELECTED FOR STUDY 


PROPERTIES AVERAGE SECTION 


SEDIMENT CHARACTERISTICS 


Min- Maxi-| Min- Min- Min- Maxi-| Mini- Maxi- Min- Maxi- 
(12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) 


reach; channel flow divided places; numerous large sand bars. Considerable winding channel; split flow; 
prevalent. Winding stream; banks vegetated sand bars visible; probably slack water areas bank 
shallow stream; sinuous; great number sand bars; split channels; vegetation negligible. Very shallow, wide 


small. 


Fairly straight channel; bed and banks all sand. Fairly straight channel; bed sand and gravel; 


13.1 (4) the bank friction comparatively negligible, and the 
active channel (even though split sand bars and possibly islands) free 
debris and vegetation. 

The computations are carried out tabular form after first determining the 
constants For practical purposes Eq. more convenient 
use than Eqs. because solutions are direct. Using the data given above 
for the problem stream, Eq. yields 


From Eq. 


The computations for three points the desired rating curve are shown 
Table procedure is: (a) assume such that Eqs. will yield velocity 
values that reasonably may expected for the given stream; determine 
and from and 19; (c) with the known value read off the value 
from the curve Fig. (d) then and (e) now 
determined and the corresponding values stage and are 
read from Fig. (f) finally, computed The relationship between 
and the desired rating curve, provided that all flow confined the 
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active channel. overbank flow occurs, then estimated allowance for 
overbank flow must made. the last column Table Manning’s 
values, computed, are shown matter interest. These, course, are 
not required for determination the relationship between and the 


TABLE DETERMINE RELATIONSHIP 
BETWEEN STAGE AND FLOW 


Stage Mann- 
(1) (la) (2a) (26) (3) (4a) (46) (6) (7) 
1.00 1.000 2.83 4.40 13.3 0.213 2.81 3.81 6.7 2,860 0.029 
2.00 1.587 4.49 2.20 20.4 0.220 3.00 5.00 8.1 6,200 0.022 
3.00 2.080 5.89 1.47 27.5 0.214 2.85 5.85 9.0 9,630 0.018 


All dimensions the foot-pound-second system. 


method described. The obvious advantage this method that the results 
not depend upon individua] judgment for selection the appropriate 
value. 
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